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O \ ABSTRACT 

We present 97 QSO candidates in 48 Galactic Bulge fields of the Optical Grav- 
itational Lensing Experiment II (OGLE-II) covering ~ 11 square degrees, which are 
O [■ selected via their variability. We extend light curves of variable objects which were de- 

tected in 3-year baseline in the OGLE-II variable star catalog to 4th year. We search 
for objects which are faint (16 < Iq < 18.5) and slowly variable during 4 years in 
this catalog by using the variogram/structure function. Finding the QSOs in such 
] stellar-crowded and high extinction fields is challenging, but should be useful for the 

astrometric reference frame. Spectroscopic follow-up observations are required to con- 
firm these QSO candidates. Follow-up observations are being prepared for four of these 
fields (BUL_SC1, 2, 32 and 45). Follow-up observations for other fields are strongly 
jjj ■ encouraged. The complete list and light curves of all 97 candidates are available in 

& ' electronic format at http://www. astro. princeton.edu~sumi/QSO-OGLEII/ 
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1 INTRODUCTION 

QSOs have been found by a variety of different methods, 
mostly using optical photometry and spectroscopy, but also 
through their radio morphology and X-ray emission proper- 
ties. 

QSO searches were usually avoided in areas where 
source crowding was a problem, for example the Galactic 
plane or the direction of nearby galaxies, such as the Magel- 
lanic Clouds. High spatial resolution of the Chandra X-ray 
Observatory limited source confusion probl ems for the X- 
ray based searches. iDobrzvcki et al.1 (l2002f) found several 
quasars behind the LMC and SMC, including four quasars 
behind the dense parts of the LMC bar by spectroscopic 
observations of the OGLE optical counterpart of Chandra 
X-ray sources. 

Another efficient method of searc hing for 
quasa rs is based on o ptical variability (lHawkins I 
1 19831: IVeron fc Hawkinsl 1995; or more recently 



Brunzcndorf & Meusingerl l2002t iRengstorf et al. I I2004D . 

Several gravitational microlensing survey groups have been 
monitoring hundreds of millions of stars toward dense 
regions such as the Galactic centr e, the Galac tic disc 
and /or the Magellan ic clouds ( EROS:lDerue Ill999t OGLE: 
Udalski et all 120031: MACHO: lAlcock et alJ l200ol : MOA: 
Bond et al]l200ll) . Everl (120021) searched QSO candidates 
towards the Magellanic clouds via variability by using 
the second phase of the OGLE exp eriment 1 (OGLE-II; 
lUjlaMcL^iubialLjL-S^Yina^ variab le star catalog 

JZebruh et al.ll200 ll) . IDobrzvcki et al. I (120031) did follow-up 
spectroscopic observations of 12 of these quasar candidates 
in the Small Magellanic Cloud fields and confirmed 5 as 
QSOs. They also found th e color information is useful to 
select QSOs. iGeha et al. 1 (120031) found about 47 QSOs 

see |http://wvw.astrou«.e du.pl/~ogle| or 

http : //bulge . princeton . edu/~ogle 



2 T. Sumi et al. 



behind the Magellanic Clouds, by selecting the candidates 
via variability in the MACHO database and follow-up 
spectroscopic observations. 

In this work, following lEver I i2002T) . we search for 
QSO candidates via varia bility in OGLE-II Galactic Bulge 
fields jlldalski et al.ll2002D . Up to now, no QSOs have been 
found in these regions. The main motivation for this study 
is that quasars may provide background sources which 
fix th e astrometric reference frame. Kinematic stu dies of 
stars JSumi. Ever fc Woznia.k^OOfllSnmi et ahlEoolh would 
greatly benefit from having absolute astrometry. These 
QSOs will also be useful for the future astrometric program, 
such as Space Interferometry Mission (SIM) and GAIA. 
However, finding QSOs in these fields is more challenging 
than in Magellanic Cloud fields because of the higher stellar 
density and higher extinction. 

In §[5] we describe the data. We select QSO candidates 
in §[3] Discussion and conclusion are given in §0] 



2 DATA 

We use the VI p hotometric maps of OGLE-II fields 
fllldalski et al] l2002L which contain VI photometry and as- 
trometry of ~ 30 million stars in the 49 Galactic Bulge fields. 
Positions of these fi elds (BUL_SC1 ~ 49) can be found in 
lUdalski et all l)2002h . We do not use BULJ3C44 in this work 
because its extinction is too high (Av > 4). The photometry 
is the mean photometry from a few hundred measurements 
in the /-band and several measurements in l/-band collected 
between 1997 and 2000. Accuracy of the zero points of the 
photometry is about 0.04 mag. A single 2048 x 8192 pixel 
frame covers an area of 0.24 x 0.95 deg. 2 with a pixel size 
of 0.417 arcsec. Details of the instrumentation setup can be 
found in Udalsk i. Kubiak fe Szvma hski (1997J). 

We also use the light curves of ob jects in the OGLE-I I 
Galactic Bulge variable star catalogue •; Wozuiak et al. 2002) 
which contain 213,965 variable objects. The variable ob- 
ject detection and the photometry have been done by us- 
ing the Difference Image analysi s (DIA) iAlard fc Lupton I 
ll998l:lAlard l200d:IWozniakl2000l) for the data collected dur- 
ing the first 3 years of the OGLE-II project. We extended 
t he light curves to the f ourth year (2000) using the method 
of lWozniak et all i2002lk however new variable object detec- 
tion is not attempted for the additional data. 



3 QSO CANDIDATES SELECTION 

We cross-refere nced 213,965 variabl e objects in the variable 
star catalogue ( Woznial^^d. 2002) with objects in OGLE- 
II photometric maps dUdalski et al1l2002l) . which are well 
calibrated, and extracted their V and / magnitudes. 205,473 
of them have counterparts in the photometric maps. We use 
these magnitudes in the following analysis. 

We chose 30,219 objects which are fainter than Iq — 16 
mag, where 7o is the J-band magnitude after extinctions 
are corrected by using the extinction map in these fields 
iSumi 1 1200411 . We show the Colour Magnitude Diagram 
(CMD) of one of these fields in Fig. |2] The detection limit 
in the variable star catalogue is about I = 19.5 mag, 
which corresponds to extinction corrected magnitudes of 



17 ~ 18.5 < Jo, depending on the level of extinction in each 
field. The nu mber density of QSO brigh ter than i = 19.1 mag 
is 9.6/deg. 2 iVanden Berk et al 112004^1 . where the Sloan i is 
almost similar to the OGLE I. So in this magnitude range 
we expect 20 ~ 30 QSOs in our 48 fields covering ~ 11 
square degrees taking the extinction into account. 

There are many artifacts in the variable star catalogue. 
The light curves of objects around bright stars are affected 
by the long wings of the Point Spread Function (PSF) of 
that bright star, and then they resemble that star's vari- 
ability. Even around a star with constant brightness, the 
wings can affect nearby faint stars due to variations of see- 
ing and the PSF. This effect causes the systematic gap in 
the brightness of faint stars around a bright star between 
in 1998 and in 1999 due to the re-aluminization of the mir- 
ror, which changed the size and/or shape of the PSF. Such 
light curves with a systematic gap are one of the major con- 
taminations in our QSO selection. So we masked the stars, 
around the bright stars including saturated stars. Saturated 
stars are not in the star catalogue. To deal with them we 
counted the number of continuous pixels, iV sa t whose ADU 
value is larger than 12,000 on the OGLE-II template image. 
Then we masked the area centered at these pixels with a 
radius of rad — 4 + 5 x V-^sat • This removed 62% of variable 
stars though only 11% of the area in total is masked. At this 
point, 11,535 l ight curves s till remain. 

Follow ing lE ver I J2 002 ) , we use t he variogram/structure 
funct ion jEver fc Gentonl Il99fll . iHuehes. Aller. fc Ailed 
1992) to find objects with increasing variability on longer 
time-scales. We measure the slope of the variogram as 
follows. For all possible pairs of measurements, differences 
of times, hij = tj — t< and square of difference of magni- 
tudes, (Ij — Ii) 2 are computed. For a given bin of log/i, the 
median of (Ij — Ii) 2 is estimated, and is noted as Var(h). 
This gives a spread of the variation of I in this time-scale. 
We bin log/i for three parts divided by 1.5, 2, 2.5, 3, then 
measure the slope, SVar = dVar(h) /dlogh. The larger S va , r , 
the larger the variability on larger time-scales. We rejected 
light curves with SVar < 0.09, which i s a sl i ghtly more 
conservative value than the 0.1 used in lEver I i2002f) . We 
search the periodicity in the light curves by using the AoV 
algorithm dSchwarzenberg-Czernv 1 ll989t) . Then we reject 
the light curves in which strong periodicity is detected with 
the AoV periodogram of larger than 50. This value was 
chosen empirically. We show two sample light curves and 
Var(h) as a function of \ogh in Fig.0 After this procedure, 
3,201 light curves remain. 

The remaining light curves, were inspected by eye. First 
we rejected apparent periodic and semi-periodic variables 
which couldn't be rejected by the previous procedure. After 
this 500 light curves are left. 

The main contamination after this selection are high 
proper motion objects and artifacts around bright Variable 
stars. We rejected them as follows. 

When an object moves from the position at the tem- 
plate image, DIA produces a pair of positive and negative 
residuals on the subtracted image. The flux of these residuals 
increases as a function of the displacement of positions of the 
objects fcver fc Wozniak Il200ll : ISoszvnski et alJl2002l). We 
cross -referenced with the proper motion catalog l|Sumi et all 
2004) to check if they have high proper motion or if there is 
a high proper motion object nearby. Even if its proper mo- 
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Figure 1. Sample light curves (first and third panel) and var- 
iogram/structure function Var(h) as a function of log/i (second 
and 4th panel, dot with solid line) for objects which passed the 
Mask around bright stars. The top two and bottom two panels 
are objects which failed (scl-2259) and passed (scl3-851) the cri- 
teria, Svar > 0.09 and AoV algorithm < 50. In the second and 
4th panels, dashed lines represent the best fit for the three points 
(log/i = 2, 2.5, 3) whose slopes Svar are written in the panels. 

tion is not so high, a bright nearby star can produces a pair 
of significant residuals and constantly increase or decrease 
the brightness of the nearby faint stars. To reject them we 
looked at the light curves of all neighboring stars in OGLE- 
II photometric maps within 20 pixels around the object, 
which were made by using DIA. Then we checked whether 
they have a companion with an opposite (increasing or de- 
creasing brightness) light curve around them. 

In this process, we also rejected the light curves which 
have neighbours with similar light curves. Usually they are 
affected by nearby bright variable stars but were not masked 
and not clearly identified in previous procedures because the 
shape of their light curves are distorted from the original 
variable due to the combination of the other artifacts or low 
S/N. The combination of the effect of the variable star and 
the high proper motion object resemble the QSO light curve. 

Finally we have 97 QSO candidates. We summarize 
our selection criteria and the number of remaining candi- 
dates after each selection criterion in Table Q We show 
20 sample light curves in Fig. [3] and list them in Table |5J 
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Figure 2. Extinction corrected Colour Magnitude Diagram of 
variable stars in BUL_SC1 (dots) and QSO candidates in all fields 
(open circles). Candidates are searched for io < 16. 

where the objects with V - I = (V - I)o = 9.999 don't 
have V-band photometry. The complete list and light curves 
of all 97 candidates are available in electronic format at 
http:/ /www. astro. princeton. edu~sumi/QSO-OGLEII/ We 
also show their position on the CMD as open circles in Fig. [3] 
As we can see in Fig. |2] many of our candidates are 
relatively bright (Iq = 16 ~ 17), which is not expected 
from a genuine quasar Luminosity Function. This means the 
contamination by irregular variable stars is higher in this 
range than expected. In this work we didn't put any cut 
by their colour, because of the following reasons. For blue 
objects, there are only f ew candidates with (V — I)o < 0.4 
JPobrzvcki et al. I j2003j) 1. The Be st ars wh ich were a major 
contamination in the previous work dEver 1 12002}) are above 
the threshold brightness Jo < 16 mag. For the redder ob- 
jects, the extinction may be underestimated for background 
objects because t he extinctions are estimated only up to the 
Galactic centre in lSumil <|2004|) . Some objects don't have V 
photometry. 
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Figure 3. 20 sample /-band light curves in 97 QSO candidates. The light curves of all 97 candidates are available in electronic format 
at |http:// www. astro. princeton. edu~sumi /qso-ogleTTTI 
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Table 2. The list of the first fO sample of 97 QSO candidates. The complete list of all 97 candidates is available in electronic format 
at http://www. astro. princeton.edu~sumi/QSO-OGLEII/ Columns arc field number, identification number as in Wozniak et al. (2001), 
ID V , and as in Udalski et al. (2002), ID, Positions on CCD in pixel, x and y, Equatorial coordinates R.A. and Dec. in equinox 2000, 
/-band magnitude and colour, I, V — I and extinction corrected magnitude and colour, Iq, (V — I)o, Variogram slope, S V ar- The objects 
with V — I = (V — I)o = 9-999 don't have V-band photometry. 



field 


ID V 


ID 


x (pixel) 


y (pixel) 


R.A. (2000) 


Dec. (2000) 


1 


V -I 


la 


(V - I) 


'S'var 


BUL.SC1 


2411 


468485 


1350.47 


4479.39 


18:02:42.91 


-29:55:02.7 


18.828 


1.195 


18.091 


0.431 


0.19 


BUL.SC2 


2329 


304696 


937.12 


3632.55 


18:04:25.85 


-28:55:46.9 


17.118 


2.614 


16.441 


1.912 


0.94 


BUL.SC2 


2423 


506014 


1282.18 


3824.72 


18:04:36.73 


-28:54:27.5 


17.416 


1.789 


16.764 


1.113 


0.68 


BUL.SC2 


3191 


337864 


947.12 


4930.96 


18:04:26.21 


-28:46:49.6 


19.155 


1.379 


18.407 


0.603 


0.15 


BUL.SC2 


3340 


136127 


369.57 


5145.97 


18:04:08.05 


-28:45:20.2 


17.081 


1.691 


16.373 


0.957 


0.13 


BUL.SC2 


4136 


165056 


138.64 


6347.10 


18:04:00.87 


-28:37:03.0 


17.087 


1.317 


16.341 


0.544 


0.43 


BUL.SC2 


4272 


164645 


385.54 


6631.55 


18:04:08.64 


-28:35:05.5 


17.289 


1.941 


16.449 


1.070 


2.31 


BUL.SC2 


4502 


775796 


1942.63 


6850.26 


18:04:57.55 


-28:33:35.6 


18.521 


1.625 


17.806 


0.883 


0.25 


BUL.SC4 


3598 


84694 


129.22 


3250.07 


17:54:07.22 


-29:49:30.2 


18.059 


1.793 


16.868 


0.558 


0.27 


BUL.SC4 


6483 


717606 


1960.77 


5624.77 


17:55:05.47 


-29:33:08.2 


18.841 


9.999 


17.697 


9.999 


0.09 



Table 1. The Number of remaining candidates after each selec- 
tion criterion. 

Criteria Candidates 



775, SC20-3616, SC32-1795 SC42-3375, SC42-3772, SC46- 
742. They have relatively constant baseline with W or V- 
shape dent in light curves. They could be new a type of 
variable star. We did not reject them because we do not 
know what they are. 

The current dataset is not optimal to find QSOs. The 
OGLE-II Galactic Bulge data are relatively shallower than 
the OGLE-II Magellanic Cloud data, which are optimized to 
find microlensing events. Variable objects were only searched 
during the first 3 years for some reason. We can improve the 
sensitivity by including the 4th year of OGLE-II and OGLE- 
III to increase the time-scale and the survey area. 

The main contamination in this work was artifacts 
around bright stars or from high prop er motion st ars rather 
than Be stars as in the previous work ySvcr 2002). Because 
the foreground stars are closer than Magellanic Clouds, most 
Be stars are above the threshold brightness. On the other 
hand, due to much higher stellar density, the artifacts are 
increased. 

The work to avoid such artifacts in the photometry on 
the subtracted images is underway. If such a new photometry 
were to be done, then the efficiency in finding QSOs would 
be increased. 
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All Variable objects 213,965 

Faint (Extinction corrected I, Iq > 16 mag) 30,219 

Mask around bright stars 11,535 

SVar > 0.09 and AoV periodogram < 50 3,201 

Eye inspection 1, cut apparent periodic variables 500 

Eye inspection 2, cut high proper motion objects 97 



4 DISCUSSION AND CONCLUSION 

We have selected 97 QSO candidates via variability from 
OGLE-II DIA variable star catalogue in the Galactic Bulge 
fields. 

Finding the QSOs in such crowded and high extinction 
fields is challenging, but should be useful for the astromet- 
ric zero-point in these fields and for the study of the fore- 
ground dust. It is especially useful for the study of the bar 
structure by using the stellar proper motions in these fields 
JSumi et al.ll2004l) . These QSOs will also be useful for the 
future astrometric program, such as Space Interferometry 
Mission (SIM) and GAIA. 

Spectroscopic follow-up observations are required to 
confirm these QSOs. The multi object and wide field spectro- 
graphs, such as VIMOS/VLT and LDSS/Magellan, are suit- 
able for this purpose. We are planning to carry out spectro- 
scopic follow-up observations in 4 of these fields (BUL_SC1, 
2, 32 and 45) by using VIMOS/VLT. Follow-up observations 
in other fields are strongly encouraged. 

As we can see in Fig. there are more candidates at 
relatively bright (Jo = 16 ~ 17) range, which is not ex- 
pected from a genuine quasar Luminosity Function. This 
means the contamination by irregular variable stars is higher 
in this range than the fainter range. So, the efficiency would 
be higher in these fainter (7o > 17) candidates. There are 
several objects which may be microlensing events, i.e. SC2- 
3191, SC2-4502, SC14-4009, SC19-4350, SC24-3042, SC30- 
3977, SC38-2298, and SC39-5030. We retained these objects 
because we do not know they are true microlensing or not. 
Also there are several objects that have light curves too simi- 
lar to each other to be QSOs, i.e. SC4-7246, SC13-332, SC13- 
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